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Abstract 

We study the conformations of a self-avoiding polymer confined to a channel by computing 
the cross-sectional distributions of the positions of its monomers. By means of Monte-Carlo 
simulations for a self-avoiding, freely-jointed chain we determine how the cross-sectional dis- 
tribution for a given monomer depends on its location in the polymer, and how strongly this 
distribution is affected by self-avoidance. To this end we analyze how the frequency of intra- 
chain collisions between monomers depends on their spatial position in the channel and on 
their location within the polymer We show that most collisions occur between closely neigh- 
boring monomers. As a consequence the collision probability depends only weakly on the 
spatial position of the monomers. Our results explain why the effect of self-avoidance on the 
monomer distributions is weaker than predicted by mean-field theory. We discuss the relevance 

of our results for studies of DNA conformations in nanofluidic channels. 
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Introduction 



The behavior of polymers in confined environments is of interest from both technological and 
fundamental perspectives. Nanofluidic channels have for example been used to stretch and vi- 
sualize single DNA molecules. ^'^ This is useful for mapping the sequence of intact long DNA 
molecules, -^"^ as well as for monitoring protein-DNA interactions,^"^ and to explore fundamen- 
tal polymer physics of DNA.^"^^ Furthermore, polymers in living systems are generally found in 
crowded and confined environments, and the confinement influences both intra- and inter-polymer 
interactions. From a polymer physics perspective, confinement affects two key parameters of 
the DNA: its orientation and its density. The orientation can be monitored by polarization sensitive 
imaging and can be used as a tool to understand the extension of the DNA. Conversely the 
local density is expected to have a direct effect on how molecules access binding sites along DNA 
where higher density gives shorter diffusion times between DNA fragments but also greater steric 
hindrance leading to lower binding rates and for example impeded overall enzymatic activity. ^^'^'^ 
Furthermore, for DNA condensation to occur, non- adjacent base pairs must approach each other. 
To obtain a detailed understanding of molecular transport and molecular interactions within con- 
fined DNA it is therefore necessary to not only know the local concentration of confined DNA, but 
also how the distribution of a section of DNA depends on its position within the molecule. Unfor- 
tunately, these monomer distributions are still difficult to determine experimentally, and analytical 
calculations and computer simulations are therefore important tools in this context. 

In the 1960s, it was realized that the monomer distributions of a confined polymer can be theo- 
retically analyzed in terms of a diffusion-annihilation equation, with absorbing boundary condi- 
tions.^^ The annihilation term represents the effect of self-avoidance, and is therefore absent for an 
ideal polymer. In this case the equation simplifies to a diffusion equation, which is easily solved 
for simple confining geometries, such as channels with circular, rectangular, or triangular ^-^ cross 
sections, or slits. Although the ideal model neglects intra-polymer interactions, it is frequently 
used to describe the monomer distribution of real polymers. ^^"^^ The reason is that calculating the 
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monomer distributions of an interacting polymer is significantly more difficult. This is because the 
annihilation term in the diffusion-annihilation equation is determined by the collision frequency, 
which in turn is related to the monomer distributions. The conventional way of approximately 
closing this system of equations is called self-consistent field theory. It assumes that the collision 
probability can be explained by a mean-field theory, neglecting intra-chain correlations. 

In this paper, we investigate the monomer distributions of a single polymer confined to a square 
channel with side length D^b, where b is the Kuhn length of the polymer. To test the effect of 
self-avoidance upon the monomer distributions and collision frequency of a polymer confined to a 
channel, we performed Monte Carlo simulations for a polymer model where the excluded volume 
of each monomer can be varied. We find that mean-field theory underestimates the probability 
of intra-chain collisions, but overestimates the extent to which this probability depends on the 
position of a monomer in the channel. This theory thus overestimates the broadening effect of self- 
avoidance. We show that the collision probability can be qualitatively understood by modeling the 
polymer as a diffusing particle with a drift velocity in the channel direction. Finally, we show 
that if the position dependent collision probability is known, the diffusion-annihilation equation 
correctly describes its impact on the monomer distributions of the self-avoiding polymer. 

Model and simulation results for the monomer distributions 

In general, the monomer distributions of a polymer confined to a channel will depend on the details 
of the polymer model, but not if the channels size (D) is much larger than the Kuhn length {b) of the 
monomer. In this limit, we can therefore restrict our analysis and simulations to the simple case of a 
self-avoiding, freely jointed polymer, consisting of N monomers of length b. Self-avoidance enters 
the model as a minimal distance a between the centers of non-neighboring monomers, yielding 
the excluded volume = Alter' jl) for each monomer. This model is illustrated in Figure 1(a). 
By changing a, the effect of self-avoidance on the monomer distributions can be modified. We 
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Figure 1 : (a) An illustration of the polymer model, showing the definition of the parameters a and 
b. (b) A snapshot (2D projection) of a configuration resulting from the Monte Carlo program, (c) 
The probability distribution of the middle monomer, Pmid(-^)> determined from simulations. The 
parameters for (b) and (c) are A'' = 800, D = 22b, a = 0.11b. 

assume that the polymer does not interact with the walls of the channel, except that the center of 
a monomer is constrained to lie within the interior of the channel. Unless otherwise stated, our 
simulations use the parameters N = 800 for the number of Kuhn length monomers in a polymer, 
and D = 22b for the channel dimensions. Figure 1(b) shows an example of a configuration, for a 
polymer with a = 0.11b. 

The cross-sectional distribution of a monomer depends on the position it occupies in the polymer. 
Figure 2 shows the distributions of the end and middle monomer of a self-avoiding chain, for two 
values of the self-avoidance parameter a. We compare the results from the Monte Carlo simulations 
to the predictions from self-consistent field theory and the ideal diffusion theory, both described 
below. While the agreement is good for small a, self-consistent theory always overestimates the 
broadening effect of the self- avoidance. 



4 



a = 0.326 a = 0.776 




x/D x/D 



Figure 2: The distribution of the end monomer (a, b) and the middle monomer (c, d) of a freely 
jointed chain ofN = 800 monomers of length b'ma channel of diameter D = 22b, for two different 
values of the self-avoidance parameter a. (a, c): a = 0.32b. (b, d): a = 0.11b. Monte Carlo 
simulations (solid line) compared to the distribution predicted by the self-consistent field theory 
(dashed line). For comparison, we also show the predicted distributions for the ideal chain (dotted 
line). 

Theoretical predictions for the monomer distributions of a con- 
fined polymer 

The starting point for understanding the monomer distributions of a confined polymer is the prop- 
agator GA?(r;ro), defined as the combined statistical weight of all polymers with one end at Tq and 
the other other end at r. If the propagator is known, the distribution of an end monomer can 
be obtained by integrating the propagator over all starting positions: pAr(r) / (iroGAf(r;ro). For 
monomer n in the interior of the polymer, the distribution can be obtained by treating the polymer 
as the concatenation of two shorter polymers, yielding p„;Af(r) oc JJroG„(r;ro) /<ir'GAf_„(r;r'). 
For this expression to be exact, the propagator corresponding to the second polymer should be 
computed in the presence of the first polymer. In the following, this complication is disregarded. 

For an ideal chain, the propagator G„(r;ro) approximately obeys a diffusion equation: 

-o'„G„(r;ro) = -b^ / evlGn{r;ro) , (1) 
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with absorbing boundary conditions at the walls. Here the operator acts on the first argument 
of G. For a polymer confined to a square channel of side length D, Eq. (1) has the well-known 
solution 

G„(r;ro) = G^ix;xo)G^{y;yo)d{z;zo), (2) 
G„ [x\xq) = 2^ ^ sm sm exp <^ \ , (3) 

k=\ V ) 

If n > (6D^) / (bn)^, Eq. (3) is greatly simplified, since the sum is dominated by the term where 
k= I. For a long chain, the distribution of a monomer at the end of the polymer and in the middle 
of the polymer therefore have the simple forms 

PendW = 2^sm — , (5) 
PmidW = ^^'^^ (6) 

These expressions are shown in Figure 2. 

It is important to consider under which conditions Eq. (1) is a reasonable simplification. For a 
self- avoiding polymer, Eq. (1) can be modified to include an annihilation term, proportional to the 
probability of collisions 

-^„G,(r;ro) = [-b^ / ev] + P„{r-rQ)\ G„(r;ro). (7) 

Pn(r;ro) is referred to as the collision probability. It is defined as the probability that a colli- 
sion would occur if one were to add a randomly oriented monomer at r to a polymer of length 
n — I, starting at tq. Eq. (7) shows that Eq. (l)-Eq. (6) are expected to work well if P„G„ ^ 
\b^/6W]Gn\ ^Pn^ 7t^b^/{6D^). 
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Figure 3: The collision probability for the end monomer of a freely jointed chain of N = 800 
monomers of length b in a channel of diameter D = 22b, for two different values of the self- 
avoidance parameter a. (a): a = 0.32b. (b): a = Q.llb. Monte Carlo simulations (solid line) 
compared to the mean field assumption Pn = ^pzPmid (dashed line). 

To bring closure to Eq. (7) in the case when this condition is not fulfilled, one must find a relation 
between G„(r,ro) and P„(r;ro). The conventional way of achieving this closure is called self- 
consistent field theory. As mentioned above, this theory rests on the mean field assumption that 
the collision probability at a given position is proportional to the local monomer density. As in the 
case of the ideal polymer, the distribution of monomer n in the interior of the chain is related to 
the propagator by p„;Ar(r) oc JJroG„(r;ro) Jdr'Gf^-„{r;r'). Since the propagator Gn approaches 
an asymptotic shape as n — )• oo, the distributions of all monomers far from the ends are expected 
to be identical. By assuming that the collision probability for a long chain is proportional to this 
universal distribution Pmid? a non-linear eigenvalue equation for the propagator can be derived and 
numerically solved (see supporting information). 




The collision probability Pn 

To test how well the self-consistent assumption for the collision probability works, we performed 
Monte Carlo simulations for a chain where overlaps between monomers are forbidden, except that 
the end monomers are allowed to overlap with any other monomer. By measuring the frequency 
of such allowed overlaps at different positions within the channel, this model allows us to measure 
the collision probability P„(r;ro) in simulations. 

In the simulation results presented here, we have only considered the cross-sectional posi- 
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tion dependence of the collision probability, which for a long chain is independent of starting 
position. We denote this quantity Pn{x,y). Figure 3 shows the marginal collision probability 
Pn{x) = J dyPn{x,y)pend{x,y) ■ The result is compared to the assumption of the mean-field theory, 
P„{x,y) = ^pzPmidi^iy)- Here p^, the number of monomers per unit channel length, is estimated 
from simulations as = A/^/ (Chain extension). Figure 3 shows that the mean-field theory un- 
derestimates the collision probability by almost one order of magnitude, yet it overestimates the 
dependence on position (x). 

That a mean-field theory for the collision probability cannot be correct is actually apparent 
from Figure 1 - this snapshot from simulations clearly shows that the average number of other 
monomers in the vicinity of a given monomer is significantly higher than one would expect from 
the average concentration p (r) . The reason is clear: since the monomers are connected, if we know 
that one monomer is located at r, closely neighboring monomers on the chain will necessarily be 
close by. Yet this correlation between neighboring monomers is neglected in the mean-field theory. 

The reason that the self-consistent field assumption concerning the relation between the col- 
lision probability and the monomer distributions does not hold can thus be understood by con- 
sidering which other monomers the end monomer most often collides with. Let us denote by 
p^^^ (r; ro) the probability that the final monomer of a chain of length A'^ collides with a monomer 
n steps removed along the chain. The total collision probability P,, is given by the sum /'„(r;ro) = 

(n) 

LnPjv (^'^o)- As before, our simulations have only measured the cross-sectional position depen- 

(n) 

dence. The n-dependence of p}^' is shown in Figure 4(a). Clearly, the collision probability is 
much higher for monomers which are closely neighboring than for distant ones. However, the 
probability of colliding with a neighbor does not vary with x in the same way as the probability 
to collide with a distant monomer does, as can be seen in Figure 4(b-d), where the collision prob- 
ability p^^\x) = J dyp^j^\x,y)pendixiy) of colliding with a monomer n steps distant from the end 
monomer is shown, for different values of a and n. In particular, for nearby monomers the collision 
probability is larger close to the walls than the self-consistent field theory predicts. 
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Figure 4: The collision probability pl^\x,y) that an end monomer with r = collides with 

a monomer n steps distant, for two different values of the self-avoidance parameter a. Blue lines: 
a = 0.32b. Green lines: a = 0.11b. Monte Carlo simulations (solid lines) compared to the collision 
probability predicted by the Eq. (10) (dashed lines), (a): p"^^ = j dxdyp"^\x,y)Psa^{x,y) . Dotted 
lines show the asymptotic scaling of Eq. (10), for large and small n. (b): p^^^ {x) / p^^\ for n = 5. 
(c): />^''(x)//>^\ averaged for n = 45 — 55. (d): /^^'■'(x)/;?^''*, averaged for n = 190 — 210. 

The collision probabilities {x,y) can be qualitatively explained. The starting point is the obser- 
vation that p^'"* must be given by 



p^^\r;ro) = / Jr'pir^(r;r';ro) ~ ^plr^(r;r;ro). 



(8) 



Here p^^' (r; r'; Tq) is the probability density that a monomer of length N, starting at Tq and ending 
at r, has monomer N — n at location r'. The integral runs over the excluded volume of the final 
monomer. pj;^\r;r';rQ) is related to the propagator of the self-avoiding chain by 



in] /■ I \ 



Giv(r;ro) 



(9) 
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Unfortunately, we do not know the propagator for the self-avoiding chain. Yet for small values of 
^ it may be a good approximation to use the ideal propagator, Eq. (3). Inserting these expressions 
for the propagator into Eq. (8) and Eq. (9), and assuming that N — n:^ Wi?' arrive at a collision 
probability which depends only on the x- and j-component of r, 

^!^^"'^) = ^VW,5u^iUJ sm-^sm-^exp| j. (10) 

Here the prefactor \J^^^ = Gn{z',z) is the return probability of the one dimensional random walk 
performed by the z-component of the ideal chain. This theory for the collision probability is com- 
pared to the results of simulations in Figure 4. Figure 4(a) shows that the theory can explain qual- 
itatively how the collision probability decreases with contour distance, although Eq. (10) underes- 

(n) 

timates the rate at which increases with separation (n). Even better, it explains quantitatively 
how the monomer collision probability depends on the position in the channel [Figure 4(b-d)], 
showing that the collision probability of nearby monomers is essentially independent of position, 
except very close to the walls [Figure 4(b)]. As the contour distance increases, the position depen- 
dence increases [Figure 4(c)]. For very large contour separations (n > 6D^/ (nb^)), the positions 
are no longer correlated, and the mean field assumption p]!j\x) °^ Pmidl-^^) ~ 2/D%vc?'{kx/D) works 
well [Figure 4(d)]. Since the probability of colliding with neighboring monomers is much higher 
than the probability of remote collisions, the total collision probability Pn = Y^hPn^ dominated 
by the terms with small values of n. Since these terms depend only weakly on position, the same 
will be true of the total collision probability. This is clearly seen in Figure 3. 

For large separations n, Eq. (10) predicts that o= n^^^^. Yet this cannot be the whole story, since 
/ n^^/'^dn diverges at infinity. Thus, Eq. (10) must break down for n above a certain threshold 
He, which we estimate by modeling the diffusion in the z-direction as a biased random walk. 
This model was developed for semi-flexible chains, but it can easily be adapted for weakly self- 
avoiding freely jointed chains. In this model, each monomer is randomly oriented, except that self- 
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avoidance causes a bias in either the positive or negative z-direction, shared by all monomers. The 
z-coordinate thus performs a biased random walk, and the probability density of any subsequence 
of n monomers starting at zo and ending at z is given by 

where /i is the bias and the variance of each step in the random walk. For broad channels, 
D ^ b, = b^/3, and the bias can be estimated by a mean field argument. The result is /i = 
[t,b^/ (6D^)] To compute the collision probability, the relevant quantity is the return probability 
density 



where 

(13) 



2o^ 2 febD 



Ai2 3 



Eq. (12) agrees with the ideal return probability y 2-^nb~ " ^ exponentially suppressed 

for large values of n. In Figure 5, the measured collision probability is compared with Eq. (10), 
modified to include the exponential cutoff. Since the length of the chain is too short to see a clear 
effect of the cutoff for the large channels we have considered so far [Figure 5(a)], we have also 
performed simulations for identical chains in a narrower channel, where according to Eq. (13) is 
smaller. For the narrower channels, it is evident that the cutoff is very well explained by the biased 
random walk model [Figure 5(b)]. 



Whereas the biased random walk model for the collision probability is qualitatively correct, it 
is not precise enough to quantitatively explain the monomer distributions of the self-avoiding 
chain. However, we can still pose the question whether the diffusion-annihilation equation cor- 
rectly describes the distribution of an end monomer, and whether the probability density of a 
middle monomer is proportional to the square of the probability density of an end monomer. To 
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Figure 5: The collision probability for the end monomer of a freely jointed chain of N = 800 
monomers, in a channel of diameter D = 22b (left) and D = 9.9b (right), for two different values of 
a. (a): a = 0.32b. (b): a = 0.11b. Monte Carlo simulations (solid lines) compared to the collision 
probability predicted by the biased random walk model (dashed lines). Blue lines: a = 0.32b. 
Green lines: a = 0.11b. Note that the collision probability is not normalized by ^, as in Figure 4. 

answer this question, we have numerically solved Eq. (7), using the value of P„{x,y) measured in 
our simulations as input. The results are seen in Figure 6, they agree very well with corresponding 
results of the simulations, even for the large value of a. As in the ideal case, the computed monomer 
distributions are somewhat too low close to the walls and too high in the center, which again can be 
explained by the fact that the diffusion approximation fails close to the wall and that both curves are 
normalized. Apart from the good agreement the most striking feature of the monomer distributions 
is how well the ideal theory works, even for the larger value of a. This is surprising, since we saw 
before that the ideal theory is expected to work well only when ^ n^b^/ (6D^). Yet Figure 3 
shows that P„ ^ 0.05 for a = 0.32b and P„ ^ 0.3 for a = 0.11b - more than one order of magnitude 
larger than n^b^/ {6D^) ^ 0.003. The resolution of this apparent paradox is that although Pn{x) is 
large, it is almost constant except very close to the walls (because most collisions are with nearby 
monomers, for which p^^\x) depends only weakly on n). Yet adding a constant term to Pn{x,y) 
only changes the normalization of the solution to Eq. (7), but does not change its shape (assuming 
P„ ^ 1 still holds). Since changing the normalization of G„ has no impact on the monomer distri- 
butions, neither has the constant term of the collision probability. Since the collision probability is 
almost constant in the interior of the channel (Figure 3), it will influence the monomer distributions 
only weakly. It is true that the sharp decrease of the collision probability close too the walls will 
cause the relative monomer density to increase significantly there, but this increase is very small in 
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Figure 6: The distribution of the end monomer (a, b) and the middle monomer (c, d) of a freely 
jointed chain ofN = 800 monomers of length b'ma channel of diameter D — 22b, for two different 
values of the self-avoidance parameter a. (a, c): a = 0.32b. (b, d): a = 0.11b. Monte Carlo sim- 
ulations (solid line) compared to the distribution predicted by the diffusion-annihilation equation 
with the measured collision probability as the sink term (dashed line). For comparison, we also 
show the predicted distributions for the ideal chain (dotted line). 

absolute terms. Thus, the effect of self-avoidance upon the monomer distributions is much smaller 
than one might expect. 



Discussion 

In this paper we have measured through Monte Carlo simulations the monomer distributions of 
a long, self-avoiding polymer confined to a square channel of size D ^ b, where b is the Kuhn 
length of the polymer. We find that the distributions are qualitatively similar to the distributions of 
an ideal chain, which can be determined from the solution of a diffusion equation with absorbing 
boundary conditions. This similarity holds even for strongly self-avoiding polymers, for which 
a priori the ideal approximation would be expected to fail. In fact, the distributions measured in 
simulations are more similar to the prediction of the ideal theory, than to the predictions of self- 
consistent field theory, which includes the effect of self-avoidance on the monomer distributions 
by means of a mean-field theory for the probability of intra-chain collisions. 

To explain this surprising finding, we have analyzed how the probability of collisions between 
two monomers depends on how close they are located along the polymer, and upon their position 
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in the channel. We have found that this collision probability can be qualitatively explained by 
calculating the return probability of a particle diffusing with a drift in the channel direction. Our 
results show that most collisions involve closely neighboring monomers, for which the collision 
probability cannot be described accurately by the mean-field theory. For these monomers, the 
collisions probability is much higher than the mean-field theory predicts, but depends only weakly 
on the position in the channel. We show that this weak dependence on position explains why 
self-avoidance has such a small impact on the monomer distributions. 

We conclude with a number of comments and questions. First, we have shown that given the 
measured collision probability P,i{x,y), the distribution of the end monomer of a freely jointed 
chain solves the diffusion-annihilation equation [Eq. (7)]. Yet to bring closure to this problem one 
must also determine how to compute the collision probability from the propagator G„, we do not 
know how to do this in a quantitative fashion. 

Second, real polymers also involve interactions between neighboring monomers. In particular, 
real polymers are usually stiffer than a freely jointed chain. In this case, the diffusion theory must 
be phrased in terms of effective monomers, each with a length corresponding to the Kuhn length. A 
case of particular interest is that of a worm-like chain with persistence length ip and effective width 
Wgff , a common model of a DNA molecule in solution. In this case, the effective monomers have a 
Kuhn length of 2ip, and the excluded volume of an effective monomer can be approximated by that 
of a cylinder, ^ = ;rweff(2^p + ^(7r + 3)£pWeff + l^^eff)-'^^ "^he results presented for a/b = 0.32 
and a/b = 0.77 would then correspond to a worm-like chain with Weff/^p = 0. 1 and Weff/^p = 0.9, 
respectively. The former value is typical for double-stranded DNA experiments performed at high 
ionic strength, ^'^^ whereas the latter corresponds to a very flexible polymer, such as single-stranded 
DNA. 32 

Third, just as the interaction between different monomers usually is more complicated than 
hard-core repulsion, so might the interaction between each monomer and the channel walls be. If 
the interaction is short-ranged and repulsive, the analysis presented here is expected to work if the 
channel diameter D is replaced by an effective diameter D^ff = dxQX^{—U {x)/kT), where U (x) 
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is the interaction energy between the walls and a monomer at position x. 

Fourth, we have seen that the self-consistent field theory cannot explain the monomer distri- 
butions of a single polymer confined to a channel (except in the limit of extremely weak self- 
avoidance, for which the ideal theory already works very well). Could it work for other geome- 
tries? In order for the monomer distributions to obey self-consistent field theory, most collisions 
must occur with distant monomers. For a polymer near a wall this is impossible, since the collision 
probability decreases faster than n^^/~, the integral of which does not converge. For a chain con- 
fined to a slit, the overlap probability of an ideal chain is marginally divergent, p°^n^^ . Therefore, 
any perturbation makes the integral over collision probabilities convergent, and long-range colli- 
sions can therefore only dominate the collision probability if the excluded volume is very small. 
For a chain that is constrained in two dimensions (a channel), we have seen that while there is a 
theoretical limit in which the self-consistent theory is valid, it is of no practical significance. 

It remains to consider three-dimensional confinement, e.g. a chain confined to a sphere or a 
cube. In this case, most collisions of a long chain can occur between distant monomers, and the 
mean-field theory might be expected to work. However, a complicating factor is that unless is 
small, a long chain will fill the available volume densely, whereas the theory is only valid in the 
dilute limit (P„ < 1). 

Fifth, what changes if instead of a single polymer, we consider a dilute blend of polymers in 
a channel? The end monomer can now also collide with monomers from other chains, and its 
distribution must change as a result. However, since self-avoiding polymers in a channel tend to 
separate, the presence of other monomers cannot influence the behavior of an interior monomer 
(unless the polymer is significantly compressed by the other polymers). Thus, the distributions of 
most monomers will be unchanged by the addition of other polymers. 

Sixth, we have seen that a self-avoiding chain in a channel can be modeled as a biased random 
walk. The strength of the bias is related to the probability of collisions, and can be estimated by a 
mean field argument. However, the mean field argument assumes that monomers are uniformly 
distributed in the channel, and should therefore underestimate the collision frequency, and hence 
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the bias. Also, since it is a mean-field theory, it does not necessarily describe the fluctuations or 

(n) 

the distribution of the end-to-end distance correctly. Therefore, the exact shape of the cutoff in p]y 
is not necessarily exponential, and the cutoff distance ric need not be related to the bias exactly as 
in Eq. (13). Another problem is that recent simulations show that whereas the biased random 
walk model works well for the monomers in the interior of the chain, monomers close to the end 
show a smaller bias. This fact should modify the probability of colliding with distant monomers, 
but the size of this effect is hard to estimate. The source of this end decay of the bias is not known, 
but might well be related to the fact that end monomers are distributed differently than interior 
monomers, as shown in Figure 6. 

Seventh, we have, throughout this paper, assumed that D ^ b, so that the diffusion approx- 
imation can be applied. It would, however, be interesting to consider the question of how the 
theory presented here fails outside this asymptotic limit. This would depend sensitively on the 
chain model, the interactions between monomers, and between the polymer and the walls of the 
channel. A particularly interesting case is that of a DNA molecule confined in a nanochannel. 
Such a molecule is usually modeled as a worm-like chain of persistence length £p ^ 50 nm, with 
screened repulsive electrostatic interactions between distant monomers, and between the molecule 
and the walls. The question of how such a molecule behaves when confined in channel sizes of the 
same order as the persistence length is of great experimental and biological interest. The resolution 
in standard fluorescence microscopy is limited by the diffraction of light, which limits the possi- 
bility to experimentally observe these distributions. However, recent experimental developments 
make it possible to directly observe DNA molecules with a resolution below the diffraction limit 
of light, and it should therefore be possible to test present and future theoretical predictions for 
the monomer distributions against experimental data. 
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Numerical solution of the self-consistent field equation 

Assuming that the density is constant along the z-direction, the self-consistent field equation ' for 
a polymer confined to a square channel is 



Xu{x^y) 
with boundary conditions 



6 V^-^^ dy^ 



+ 3^ ) +^Pzu^{^.y) 



u{x,y), (1) 



u{x,0) = u{x,D) = u{0,y) = u{D,y) = 0. (2) 
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Here, the eigenvalue A is a free parameter, the parameters b and ^ are given by the polymer model, 
and is estimated from simulations as = N/ (Chain extension). Since this eigenvalue equation 
is non-linear, it is important to note that the eigenfunction u{x^y) must be normalized, 

Jj'^dxdy[u{x,y)f = 1. (3) 



We have solved this equation numerically by an iterative procedure. Starting with uo{x,y) = 
{2 / D) s'm{nx / D) sm{Ky / D) (the ground state of the diffusion equation), we updated the solution 
according to 



1 



[{\-a)u„ + au'„] 



(4) 



where Zn+i is a normalization constant. is the normalized solution to the (linear) eigenvalue 
equation 



Xu„{x,y) 



+ 



+ <^p2 [un{x,y)X 



u'„{x,y). 



(5) 



6 \ dx^ dy^ 

The convergence constant a is necessary for the algorithm to converge to a solution. For the 
parameters we have considered, a = 0.5 is enough to ensure convergence. 

Eq. (5) was solved by discretizing the region <x,y < D/2 into a grid of 100 x 100 squares and 
representing the function u',^{x,y) as a vector of length 100 x 100 = 10000. Treating this region 
suffices, since the other three quarters are determined by symmetry. In this representation, the 
operator 



+ ■ 



+ ^p.[unix,y)y 



(6) 



6 \ dx^ dy^ 

with appropriate boundary conditions, is a sparse matrix of size 10000 x 10000. The eigenvector 
corresponding to its smallest eigenvalue gives the solution m',(jc, j). This eigenvector was deter- 
mined by using the subroutine eigsh of Numerical Python (NumPy).^ 
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Details of the Monte Carlo program 

Our Monte Carlo simulations implement the Metropolis algorithm, with crankshaft trial updates. ^ 
In each update step, n adjacent molecules / = nmin, Hmin, . . . , n^ax = «min + n — 1 are selected, where 
1 < n < 16. If the n monomers are all in the interior of the chain, they are rotated by a random 
angle 6 G [0, In) around the axis connecting monomers Hmin — 1 and /imax + 1 . If instead Hmin = 
1 (or Hmax = N) they are rotated by a random angle E [0, 2n) around a random axis passing 
through Hmax + 1 (or Hmin — 1). Note that these trial updates conserve the distance between adjacent 
monomers. For these simulations, the chain is initialized such that the distance between the center 
of adjacent monomers is b. In order for the chain to equilibrate, 2 • 10^*^ trial updates are suggested 
before starting the measurement process. The measurement process runs for 10^^ trial updates, 
with measurements taken once every 10^ steps. 
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